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We show that binary mixtures of Bose-condensates of alkali atoms have a 
great variety of ground state and vortex structures which can be accessed ex- 
perimentally by varying the particle numbers of different alkalis. We have con- 
structed a simple algorithm to determine the density profiles of these states, 

00 ' and have worked out their phase diagrams within Thomas-Fermi approxima- 

O ' 

tion. Many structures of the alkali binary contain a coexisting region, which 

On 



is the analog of the long sought 3 He- 4 He interpenetrating superfluids in ultra- 



& ' low temperature physics. 

The search of Bose condensate in alkali atoms Q || Q has a deep root in ultra-low 
temperature physics. Since the discovery of superfluid 3 He, the searches of the next elemental 
superfluid have been focusing on spin polarized hydrogen and 3 He- 4 He mixture. The former 



^ ■ promises another Bose superfluid besides the only known example of 4 He, the latter, the first 



example of interpenetrating superfluids. The recent discoveries of alkali Bose condensates 
@ have in essence achieved the goal of the superfluid hydrogen search. Since there 
are no intrinsic difficulties in loading more than one alkali element and have them cooled in 
the same trap, it appears highly promising that interpenetrating superfluids may be realized 
for the first time within the same experimental setting. 

In this paper, we shall discuss binary mixtures of alkali condensates. Such mixtures 
may consist of different alkalis such as 87 Rb- 23 Na, or different isotopes such as 87 Rb- 85 Rb, 
or different hyperfine states of the same alkali such as the (F = 2, Mp = 2) and (F = 
l,Mp = 1) states of 87 Rb. We shall denote the two different alkalis as 1 and 2, and their 
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particle numbers as JVi and N 2 . Unlike single component systems which are characterized 
by a single scattering length, alkali binaries are characterized by three scattering lengths 
ai,a 2 , and 012; representing interactions between like and unlike alkalis. (While scattering 
lengths of like alkalis are known at present, those of unlike alkalis are not). As we shall see, 
this moderate increase in energy scales leads to a proliferation of ground state and vortex 
structures, which we shall illustrate for the case ai,a 2 > 0. This case is chosen because it 
has the greatest structural diversity, and can be analysed by simple analytic methods. The 
results for this case will also be useful in understanding the qualitative features of other 
(negative scattering length) cases. 

Our key results are : (a) a simple algorithm for determining the density profiles of the 
binary mixtures, (b) the phase diagrams of the vortex free ground states [denoted as (vO)], 
and the vortex states in either 1 or 2 [denoted as (vl) and (v2) respectively]. We limit 
our discussions to these two types of vortices because they are the states the system will 
first fall into as the confining potential is rotated. Our algorithm, however, can be applied 
to arbitrary number of vortices in 1 and 2. As we shall see, the structure of the mixture 
depends on two parameters a and (3, which are proportional to the strength of a 12 and the 
relative strength between a± and a 2 respectively. These parameters determine whether alkali 
2 when added to an existing cloud of 1 will stay at its exterior or interior. Many structures 
of the mixture contain a coexisting region which is especially large when Ni ~ N 2 . This is 
the analog of the long sought 3 He- 4 He superfluid mixture in ultra-low temperature physics. 
Even greater structural diversity is found when vortices are inserted in one of the alkalis. 
When the alkali 1 and 2 are mutually repulsive, a 12 > 0, the vortex free alkali (say, 2) always 
enter the vortex core of 1, giving rise to a variety of "vortex donut" structures in 1. When 
a 12 < 0, one has an unusual "concentric donut" structure where a vortex free "donut" of 2 
embedded in a "vortex donut" of 1. 

Our results are obtained by minimizing the Gross-Piteavskii energy E(fy 1 , ty 2 ) = T + V 
subject to the constraint of constant particle numbers, i.e. by the condition 5K = 0, 
K = .E^i, ^2) — A^i Ni —fi 2 N 2 , where //j) are the order parameter and chemical potential 
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of the i-th alkali, % = 1,2. The potential energy V and kinetic energy T are 

Mi 1 i y 



/f«(x)i«,i" + ^i»,i'U/^i« l H« J |" d) 



i=l,2 



T = / ( E - + Ci*iv*i • *;w 2 + c 2 ^v^i • # 2 w* + c.c. ) . (2) 



v i=l,2 

C/j is the potential of the magnetic trap for the i-th alkali, generally of the form [|]] [Zj(x) = 
3 ' 2 f ( r2 + A 2 ^ 2 )= \MiUjf{r 2 + A 2 z 2 ), where g»j and M* are the g-factor and mass of the i-th 
alkali, /i_b is the Bohr magneton, B a is the magnetic field at the center of the trap, L is the 
length scale of the variations of the magnetic field, and A is the trap anisotropy. £i and C, 2 
are complex coefficients which is caused by back flow effects between different alkalis. They 
are expected to be small in the dilute limit. 

Our calculations are performed within the Thomas- Fermi approximation (TEA), which 
is a good approximation in the limit of large number of particles, as pointed out by 
Baym and Pethick ||. In TFA, one ignores all V|^t| terms in T. For vortex free 
structures (vO), this amounts to setting T = 0. For (vl) vortices, = |\l/i|e^ and 
^2 = 1^21) , TFA amounts to retaining only the centrifugal term |\l/i| 2 /r 2 in T. Be- 
cause of the absence of gradient terms V|^i|, the densities obtained from TFA usually 
consist of discontinuities in radial curvature. Such changes in curvature in fact take 
place smoothly over the distance of correlation length, which decreases with increasing 
particle numbers f|. It is convenient to write K and Ni in dimensionless form by the 
following rescaling: (i) z — > z/X, (ii) x — > a^x, where h 2 jM\a% = gi/i#I? /L 2 , (iii) 
|^i| 2 = Aift, where A 1 = M^/iBnhai), A 2 /A x = (a 1 M 2 /a 2 M 1 ) 1 / 2 , (iv) p x = 5^1^ 
and /i 2 = \%Ui{A 2 /A x )v 2 , (v) K = {fiui/ AX){a T / {knaijK,. (vi) N\ = (87rA) _1 (ar/ai)ni, 



Ui= d 3 xpi (3) 



N 2 = (87rA) 1 (a T /ai)^aiM 2 /a 2 Min 2 . We then have 

{ r 1 

K = J d 3 x -Cipi - C2P2 + ~ (pl + p\ + 2a PU 

ci(x) = v x - (r 2 + ^ + z 2 ^j c 2 (x) = v 2 - (3(r 2 + z 2 ), (4) 



where x = (r, z) = (x,y,z) in eq.(|]), and p = and 1 for (vO) and (vl) state respectively. 
a and (3 are defined as 

«12 a 92 la 1 M 2 , . 

a = , P=—\—n-> 5 

reflecting the interaction between unlike alkalis, and the relative strength of the self interac- 
tion between like alkalis. For (v2) vortices, the energy can be rescaled to the same form as 
eq.(||) with the same a and (3 in eq.@, and with (ci=r 2 + z 2 — ux, c 2 =(3{r~ 2 + r 2 + z 2 ) — v%) 
||. We shall not consider the case (a < —1), which does not have bulk stability. In our 
subsequent illustrations, we use the example of 87 Rb- 23 Na mixture, (our alkali 1 and 2), 
which has identical g factors. We thus have ax ~ 100 A |L], a 2 ~ 49 A ||, and (3 = 0.73. 
Since a\ 2 is not known for this mixture, we consider the entire range of a > —1. Taking 
A = 1 and a? = 4 x 10~ 4 cm, we have Nx = 16.7ni, N 2 = 12.3n 2 - 

Our goal is to find the condensate structure (i.e. px, p 2 ) as a function of particle number 
nx,n 2 - This is done by : (I) Minimizing K, to find the equilibrium densities pi for given 
chemical potential z/j, (II) substituting these densities into eq.(|3]) to obtain the relation 
rij = riiiyx, v 2 \ (III) inverting this relation to obtain z/j = z^(ni, n 2 ) and hence the evolution 
of pi (through their dependence on z/j) as a function of rij. Although much of our labor went 
into (II) and (III), they are straightforward (though lengthy) calculations once the densities 
profiles are determined by the simpler but subtler step (I), which we now explain. 

Let [0], [1], [2], [12] denote the vacuum, the single phase of 1, 2, and the coexisting phase 
of 1 and 2 respectively. Their densities are given by the stationary conditions of /C, 

[12]: Pi = ^?>0, p 2 = ^^>0; [0]: Pl = p 2 = 0; (6) 

1 — a x — a 

[1] : Px = cx > 0, p 2 = 0, [2] : p 2 = c 2 > 0, Pl = 0. (7) 

The distribution of these "phases" in Cx-c 2 plane will be referred to the distribution plot, 
which depends only on a. The distribution plots for (0 < a < 1) and (— 1 < a < 0) 
are shown in fig. (1.1) and fig. (1.3). [The distribution plot for (a > 1) will not be shown 
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as it is given by fig. (1.1) with the coexisting region [12] collapsed into the line c\ = c 2 .] 
The boundaries separating [12] from [2] and [1] are denoted as 1 Q and 2 Q . They are the 
surfaces of vanishing p\ and p 2 , described by equations (1 Q : C\ = ctc 2 ) and (2 Q : c 2 = aci) 
respectively. The boundaries separating [0] from [1] and [2] are denoted as 1° and 2°. 
They are described by equations (1° : C\ = 0) and (2° : c 2 = 0). ^From eq.(||), one can 
see that a path in real space will have an "image" path in Ci-c 2 space. For example, the 
image of a radial path on the horizontal plane with height z is given by T(z) : C\ — v\(z) = 
/3 _1 [c 2 — 1*2(2)] + pP[c2 — V2{z)]~ l , where v\{z) = v\ — z 2 , vi{z) = 1^2 — /3z 2 . For vortex free 
state (p = 0), T(z) is a straight line with slope f3 emerging from the point (ui(z), v 2 (z)), 
which is shown as a dashed line emerging from a circle in fig. (1.1). For the (vl) vortices, 
{r(z)} is a family of curves shown in fig. (1.2), where Z2 > Z\ > in this figure. The arrows 
on these paths indicate the direction of increasing r. As r varies from to 00, T (z) intersects 
the phase boundaries in a specific order. ^From the definition of Vi(z), it is easy to see that 
as |2:| increases, T(z) slides down rigidly along the straight line with slope f3. 

The above considerations suggest a simple algorithm for determining the structure of 
the alkali binary : (i) For given z/!,z/ 2 , draw the image paths {r(,z)} on the distribution 
plot. T(z) will intersect a set of phase boundaries in a specific order, (ii) Construct the 
set of boundary surfaces intersected by {r(z)}. (iii) Eliminate all portions of the boundary 
surfaces that are inconsistent with the order of intersections generated in (i). The remaining 
surfaces are the physical boundaries in the mixture. The densities of various phases bounded 
by these surfaces are given by eq.(H) and (0). 

To illustrate this algorithm, consider the vortex free mixture in fig. (1.1). The image 
paths of this mixture intersect phase boundaries 2 Q , 1 Q and 2°, corresponding to a set of 
concentric spherical surfaces in real space separating the single component and coexisting 
regions as shown in fig.(1.4). Fig.(1.4) is the structure of the mixture. For the (vl) vortex in 
fig. (1.2), the family {r(z)} intersects all four phase boundaries (1 Q , 2 Q , 1°, 2°) in the order 
(1 2 1°), or 2° alone. Displaying all four boundary surfaces as in fig. (1.5) and eliminating all 
sections of these surfaces inconsistent with the order of intersection, we obtain the physical 
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boundaries shown in fig. (1.6). The densities of the phases [1], [12], and [2] enclosed by these 
surfaces are given in eq.(||) and (0). It is clear that this algorithm applies to mixtures with 
arbitrary number of vortices in 1 and 2, which have different image trajectories. 

Having determined the density profiles for given chemical potential (1/1,1/2), we have 
followed the aforementioned Steps (II) and (III) to construct the phase diagrams for both 
ground states [(vO)] and vortex states [(vl) and (v2)] over the entire range of (a,j3). For 
reasons to be explained shortly, we shall focus on the parameter range [A: < a < 1, /3 < a], 
[B: < a < 1, a < (3 < 1], and [C: — 1 < a < 0, (3 < 1]. For brevity, we shall denote the 
(vO) state of A as (vO)-A, the (v2) state of B as (v2)-B, etc. As it turns out, for ground 
state phase diagrams, it is sufficient to discuss those of (vO)-A and (vO)-B. For vortex 
phase diagrams, it is sufficient to discuss those of (vl)-A, (v2)-A, (vl)-B, and (vl)-C. 
The phase diagrams of all other cases can be obtained from these either by interchanging 1 
and 2 or by collapsing the coexistence region down to a line . In all cases we have studied, 
the condensate structure evolves continuously over the entire ni-n 2 space, even though the 
phase diagrams may contain different structural regimes. 

The typical phase diagrams of (vO)-A and (vO)-B are shown in fig. (2.1) and (2.2), with 
(a = 0.9, (3 = 0.73) and (a = 0.6, (3 = 0.73) respectively. The density profiles for mixtures 
marked by letters "a" to u g" in fig. (2.1) and (2.2) are shown in fig. (2.3) to (2.9). The line 
between b and c in fig. (2.1) divides the configurations where 2 is absent or present at the 
origin. The main difference between A and B is that the former (f3 < a) represents a regime 
where repulsion between unlike alkalis dominates. As a result, when a small amount of 2 is 
added to an existing condensate of 1, it stays at its surface in the case of A (see fig. (2.3) and 
(2.4)), whereas it enters directly into the center of 1 in the case of B (see fig. (2. 7)). One can 
see from fig. (2.5) and (2.8) that when 1 and 2 have similar number of particles, the coexisting 
region of 1 and 2 occupies a substantial portion of the mixture. This large coexisting region 
is the analog of the long sought interpenetrating superfluid phase in 3 He- 4 He mixture. 

Next, we turn to the vortex states. The presence of a vortex (say, in alkali 1) will turn 
the spherical ground state condensate into a vortex "donut" . The donut hole is the vortex 
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core, around which supercurrent circulates. For repulsive interactions, (a > 0), the general 
behavior of the (vortex free) alkali 2 is to fill up the vortex core and to penetrate into the 
interior of vortex 1. The typical phase diagram of (vl)-A, (v2)-A, and (vl)-B are shown 
in fig. (3.1), (3.5), and (3.7) respectively, with values of (a = 0.9, (3 = 0.73) for A and 
(a = 0.6,(3 = 0.73) for B. The scales of n% and n 2 are chosen so that a condensate with 
N\ = N 2 = 10 4 will appear in the middle of the diagram. There are in fact many more 
structural regimes near the n\ and n 2 axes which correspond to various stages of filling up 
the vortex core. However, they only exist close to the n\ and n 2 axes (i.e. small number of 
particles in either 1 or 2) and are not visible on this scale ||. 

The boundary surfaces for the mixtures marked as "a" to "/" in fig. (3.1), (3.5) and (3.7) 
are shown in the rest of figure 3 with self explanatory labelling. These surfaces are presented 
only in a quadrant of r-z plane, as they are cylindrical symmetric about z and have mirror 
symmetry about the xy plane. The boundaries 1 Q and 2 Q are represented by solid lines, 
while (1° and 2°) are represented by dashed lines. 1 and 2 coexist in the region between the 
dashed lines. The curve in fig. (3.1) between a and b divides regions where a ring of single 
phase [1] is present or absent in the interior of the "vortex donut" . The curve between b and 
c divides the regions where the "vortex donut" resides entirely inside or extends beyond the 
spherical cloud of 2. The phase diagram of (v2)-A is very simple as all mixture states have 
the same qualitative structure shown in fig. (3.6). The phase diagram of (vl)-B (fig. (3.7)) 
is essentially that of (vl)-A (fig. (3.1)) with the region occupied by "6" is collapsed into a 
line. The absence of this region is due to the weakness of mutual repulsion between different 
alkalis (i.e. small a), so that alkali 2 always penetrate completely the vortex donut of 1. 

Finally, we consider the (vl)-C vortices, i.e. the vortex state of a mixture with a\ 2 < 0, 
(hence a < 0). The typical phase diagram is shown in fig. (4.1), (where we have chosen 
a = —0.9,(3 = 0.73). It is qualitatively the same for all values of (3. The boundary surfaces 
of the mixtures marked as a, b, and c are shown in fig. (4.2) to (4.4). The region in fig. (4.1) 
containing "a" describes the unusual structure of a vortex free donut of 2 is embedded in 
the vortex donut of 1. When crossing the boundary separating a and b in fig. (4.1), the hole 
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of "donut" 2 shrinks to zero, filling up the vortex core of f. As one crosses the boundary 
separating b and c in fig. (4.1), the "vortex donut" 1 is completely devoured by 2. 

^From the above discussions, we see that binary mixtures of alkali Bose condensates 
possess a great variety of ground state and vortex structures, which can be accessed by 
varying the number of particles of each alkali. This array of structures allows one the go 
continuously from a regime of interpenetrating superfluids to one with separated phases. The 
possibility of scanning through this continuum offers great opportunities to study coupled 
macroscopic quantum phenomena and interactions of elementary excitations in distinct Bose 
fluids. Realizations of alkali condensate mixtures will certainly deepen our understanding 
of interacting Bose fluids, and widen our horizon on superfluid phenomenon. 

TLH would like to thank Greg Lafyatis for discussions. This work is supported in part 
by NSF Grant No. DMR-9406936. 
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Figure 1: Fig. (1.1) and (1.2) show the typical distribution plots for < a < 1. The "image 
trajectories" T(z) (dashed lines) of the ground state (vO) and the (vl) vortex of a mixture 
with chemical potential (ui, v 2 ) are shown as dashed lines in fig. (1.1) and (1.2). In fig. (1.2), 
z 2 > z 1 > 0. Fig.(1.3) shows the distribution plot for -1 < a < 0. Fig.(1.4) and fig. (1.5) 
to (1.6) show the real space boundary surfaces for (vO) and (vl) cases respectively. The 
boundary surfaces infig.(1.4) are given by [2° : r 2 +z 2 = v 2 / f3\, [1 Q : r 2 +z 2 = (ui — au 2 ) / (1 — 
a/3)], [2 Q : r 2 + z 2 = (av\ — u 2 )/(a — (3)}. The boundary surfaces in fig. (1.5) are given by 
[1° : r 2 + r~ 2 + z 2 = [2° : r 2 + z 2 = v 2 / (3\, [1 Q : (1 - a(3)(r 2 + z 2 ) + r~ 2 = v x - av 2 \, 
[2 Q : (a — j3)(r 2 + z 2 ) + a/r 2 = av\ — v 2 \. Fig. (1.6) gives the structure of the mixture in the 
r-z plane. 

Figure 2: Typical ground state phase diagrams for (vO)-A and (vO)-B. We have chosen 
{a = 0.9, (3 = 0.73) for A, {a = 0.6,/? = 0.73) for B. The dotted line in fig. (2.1) represents 
the case of equal particle numbers in both condensates, (iV\ = N 2 ). As discussed in the 
text, for parameters appropriate for 87 Rb and 23 Na, and for current magnetic trap, we have 
Ni = l6.7ni,N 2 = 12.3n 2 . The circle in fig. (2.1) indicate a mixture with Ni = N 2 = 10 4 . 
For the densities shown in fig. (2.3) to (2.9), the chemical potential v for the outer cloud is 
of the order of 10 to 20, giving rise to a mixture size of about 3 to 4 in dimensionless units, 
which is about 2 x 10" 3 cm for a cloud of N x ~ 10 4 , A^ 2 ~ 10 4 . 

Figure 3: Typical phase diagrams for (vl) and (v2) vortices for parameter range [A : < 
a < 1, f3 < a], and [B : < a < 1, a < f3\. The values of (a, (3) used here are identical to 
those used in figure 2. Figures with letter labels "a", "6", etc. at the upper right corners 
show the boundary surfaces of the mixtures indicated by the same letter in the n\-n 2 phase 
diagrams. 

Figure 4: Fig. (4.1) shows the typical phase diagram for the (vl) vortices for case [C: 
a < 0], with (a = —0.9,(3 = 0.73). The notations here are identical to those in figure 3. 
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